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ABSTRACT 

We have coupled a chemical model with two dynamical models of collapsing low mass 
star-forming cores to predict abundances across the core of the commonly used infall 
tracers, CS and HCO"*", at various stages of the collapse. The models investigated are 
a new ambipolar diffusion model and the 'inside-out' collapse model. We have then 
used these results as an input to a radiative transfer model to predict the line profiles 
of several transitions of these molecules. For the inside-out collapse model, we predict 
significant molecular depletion due to freeze-out in the core centre, which prevents the 
formation of the blue asymmetry (believed to be the 'signature' of infall) in the line 
profiles. Molecular depletion also occurs in the ambipolar diffusion model during the 
late stages of collapse, but the line profiles still exhibit a strong blue asymmetry due to 
extended infall. For the inside-out collapse model to exhibit the blue asymmetry it is 
necessary to impose a negative kinetic temperature gradient on the core and suppress 
freeze-out. Since freeze-out is observed in several class protostars which are thought 
to be collapsing, this presents a major inconsistency in the inside-out collapse model 
of star formation. 
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1 INTRODUCTION 

The current method for detecting infall in pre-stellar cores 
and young stellar objects (YSOs) involves observing an op- 
tically thick line, such as CS J=2— >1 and 3— >2, along with 
an optically thin line, such as N2H"'" J=l— >0 (both of which 
are tracers of dense gas) towards the selected core (Lee et 
al. 2004). If the optically thick lines exhibit a double peaked 
profile, with the blue-shifted peak stronger than the red- 
shifted peak (i.e. a blue asymmetry), and the optically thin 
line is single peaked with a central velocity that lies inbe- 
tween the two peaks of the optically thick transition (to 
verify that the double peak does not arise from two sepa- 
rate dense clouds along the line of sight), then it is a strong, 
but not unambiguous, indicator that the core is collapsing. 
The blue asymmetry is not unique to a particular col- 
lapse model (Leung & Brown 1977), although many as- 
tronomers interpret their observations of both spectral line 
and continuum emission by comparison with the inside-out 
collapse model of Shu (1977) (e.g. Zhou et al 1993; Motte 
& Andre 2001). The inside-out collapse model is also widely 
used in theoretical studies of collapsing cores (e.g. Tsamis 
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et al. 2008; Rawlings & Yates 2001; Rawhngs et al. 1992; 
Evans et al. 1994; Weidenschilling et al. 1994). 

However, there are several theoretical arguments and 
observational studies which indicate that the inside-out col- 
lapse model may be inappropriate. For example, IBelloche] 
|et al.| ( |2002[ ) found that the class protostar TRAM 04191 
shows evidence of extended inward motions, which are in- 
consistent with the model of Shu (1977), which predicts that 
the outer parts of the core should be static. 



The inside-out collapse model of Shu ( 1977 1 describes 



the protostellar phase, and begins from a singular isother- 
mal sphere. Preceding this is the pre-protostellar phase (of- 
ten referred to as the pre-stellar phase for breivity; |Ward-| 



Thompson et al. (1994l). It has been suggested that during 



this time a diffuse core evolves quasi-statically into a more 
centrally condensed pre-stellar core via the process of am- 
bipolar diffusion. In this paper, we model the line profiles of 
various transitions of the commonly used infall tracers CS 
and HCO^ produced by two different models of collapsing 
cores: (i) a new ambipolar diffusion model (Stace & Rawlings 
2010, in preparation; hereafter SRIO) and (ii) the standard 
inside-out collapse model (Shu 1977). The ambipolar diffu- 
sion model covers the pre-stellar phase, but by ~ 10® yr it 
extends into conditions similar to the protostellar regime, 
reaching a density of nn ~ 2 x 10^ cm~^ at a distance of 
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10^® cm (~ 700AU) from the core centre, which is actually 
higher than that estimated for the protostellar core IRAM 
04191 at such a distance ( Belloche et al.|2002 | At this time, 
one of the main differences between these two models is that 
the inside-out collapse model has a static envelope, whereas 
in the ambipolar diffusion model there are extended low- 
velocity inward motions. 

An important aspect of this work is that we couple the 
dynamical models to a large chemical network, in order to 
predict more accurate abundance profiles of each molecule. 
These abundance profiles are then used as an input to a ra- 
diative transfer model, along with the density and velocity 
profiles appropriate for each collapse model. We focus on 
how the chemical evolution of the cores affects the resulting 
line profiles shapes, and we find that the chemical distribu- 
tion can be an important factor in determining whether or 
not the blue- asymmetry is seen. 

The paper is organised as follows. In Section [2] we de- 
scribe dynamical collapse models studied in the paper, and 
in Section [3] we describe the chemical model, the results of 
which are given in Section [4] In Section [S] we describe the 
radiative transfer model used and how the line profiles were 
generated, and in Section [6] we present the resulting line 
profiles, and also discuss their sensitivity to the tempera- 
ture profile of the cores. In Section [T] we present a discus- 
sion of the results, considering their sensitivity to the con- 
ditions chosen in the chemical model. The conclusions are 
summarised in Section [S] 



2 DYNAMICAL MODELS 

2.1 Ambipolar diffusion 

Ambipolar diffusion controlled collapse occurs in magnet- 
ically subcritical cores, where the magnetic field is strong 
enough to support the core against gravitational collapse. 
It is estimated that a core is subcritical if the core mass to 
magnetic flux ratio, M/(j), is less than ~ O.ISG"^" (in cgs 
units), where G is the gravitational constant (Mouschovias 
& Spitzer 1976). In such cores, the neutral species do not 
directly feel the effect of the magnetic field, and therefore 
drift towards the centre of core, being impeded by the drag 
force exerted by the ions which are tied to the magnetic field. 
Eventually, in the centre of the core, enough mass builds up 
so that the core becomes supercritical and begins to collapse 
on a dynamical timescale. 

We have taken preliminary results from the dynamical 
model from SRIO to couple with our chemical model. The 
model will be described in detail in SRIO, but the following 
is a brief summary: 

The dynamical model follows the spherically symmet- 
ric collapse of a core composed of 100 concentric shells of 
gas and dust. The core is embedded in a cylindrical mag- 
netic field, which evolves throughout the collapse, although 
its geometry always remains cylindrical. There is no mass 
transfer between the shells. 

At each time step, the code calculates the total force 
exerted on each shell of gas, taking into account the grav- 
itational force, the gas pressure and the drag force exerted 
on the neutral species by the ions and charged dust grains. 
These drag forces depend on the number density of charged 



particles, and on the drift velocity between ions and neu- 
trals, which is related to the magnetic field strength. To 
calculate the magnetic field at each step, it is assumed that 
the magnetic field is frozen into the plasma. 

The code also calculates the visual extinction for each 
shell, which increases during the collapse due to the increase 
in the density of dust grains. A relationship between visual 
extinction and temperature is assumed, so the radial tem- 
perature profile also evolves throughout the collapse, which 
affects the gas pressure. The temperatures in the core ranges 
from 6.5 K to 13.5 K, although since in our chemical model 
it is not yet possible to include a time dependent chemistry, 
we have assumed a constant gas temperature of 10 K. The 
chemical abundances are not significantly affected by this 
simplification; running the chemical model at temperatures 
of 6.5 K and 13.5 K results in differences in the CS and 
HCO+ abundances of < 30%. 

In the model of SRIO, account is taken of how the chang- 
ing chemical composition of the core affects the collapse dy- 
namics; chemical reactions and freeze-out of ions onto dust 
grains affect the overall ionisation fraction of the gas, which 
therefore alters the drag force experienced by the neutrals. 
The rates of the chemical reactions depend on the density, 
visual extinction and temperature calculated by the dynam- 
ical equations. SRIO have coupled the chemical and dynam- 
ical equations to build a self-consistent model of ambipolar 
diffusion controlled collapse. 

We have used the results for a core of mass 6.80 Mq, 
radius 8.45 x 10^^ cm and magnetic field strength 4 
/iG, leading to a mass to flux ratio of M/{B-kF&) ~ 
1500 g G~^ cm~'^. This is ~ 3 times the critical value cal- 



culated by Mouschovias & Spitzer (19761, consistent with 
observations of molecular clouds ( [Crutcher et al.||1999[ ). 
Although this is a 'supercritical' core according to the 
Mouschovias &: Spitzer] (1976) criterium, the magnetic field 



is still strong enough to significantly retard the collapse. The 
initial density distribution of the core is uniform, with wh = 
2.8 X 10'^ cm~'^. The data we have taken from this model are: 
(i) radius as a function of time, and (ii) density and velocity 
as a function of time and radius, for the set of 100 infalling 
shells making up a collapsing core. We run chemical models 
for six of these shells, spanning the entire core radius. We de- 
rive the chemical evolution of the shells with starting radii of 
8.45 X 10^^ cm, 7.37x 10^'' cm, 6.77x lO" cm, 5.41 x 10^'^ cm, 
3.93 x 10^'^ cm and 1.82 x 10^^ cm respectively. 

The physical evolution of the test shells is shown in 
Figure IT] The core contracts with time (see the evolution of 
the outer shell, plotted with the black solid line in Figure [1]), 
as the density of the inner shells increases. The inner shell 
(grey dashed line), collapses relatively quickly into the core, 
and before 10^ yr it reaches its final density of ~ 10^ cm"'' 
where the collapse halts for this shell. This is because for 
such high and rapidly changing densities, the model of SRIO 
cannot converge on a solution. 

Since the dynamics of SRlO's model depend upon the 
chemistry, particularly the ionisation fraction of the gas, it 
is possible that a different set of chemical reactions and 
initial conditions could be inconsistent with the dynamics 
predicted by SRlO's models. The set of reactions used in 
the chemical model of SRIO was taken from the UMIST 
RATE95 database (Millar et al. 1997), whereas in our model 
we have used a more up-to-date reaction set from the 
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Table 1. Table showing the initial radius (i?o) and density (nno) 
for each collapsing shell for which we derive the chemical evolu- 
tion, for the ainbipolar diffusion and inside-out collapse models. 



Ambipolar Diffusion 


Inside-out 


Total mass 6.80 Mq 
Outer radius 8.45 X 10^'^ cm 


Total mass 0.96 Mq 
Outer radius 1.60 X 10^'^ cm 


i?0 (cm) riHO (cm-3) 


Ro (cm) nno (cm'^) 



1.82 X 10^^ 
3.93 X 10^^ 
5.41 X 10^^ 
6.77 X 10^^ 
7.73 X IQi'^ 
8.45 X IQi'^ 
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Figure 1. Physical evolution (radial position, velocity and den- 
sity) of the test shells in the ambipolar diffusion model. The black 
solid line, black dotted line, black dashed line, grey solid line, grey 
dotted line and grey dashed line are for the shells beginning at 
Ro = 8.45 X 10^^, 7.73 x 10^'', 6.77 x lO^'^, 5.41 x 10^^, 3.93 x lO^'^ 
and 1.82 X lO^'^ cm respectively. 



RATE06 database (Woodall et al. 2007) (see Section |3|. 
Furthermore, SRIO calculated the visual extinction self con- 
sistently with grain distribution in the core, whereas for sim- 
plicity we assumed a constant visual extinction, which will 
have an effect on the rates of photoreactions, particularly at 
the edge of the core. 

Our chemical model predicts an ionisation fraction of 
~ 10~® throughout the core at early times (~ 10^ yrs), 
which is approximately an order of magnitude higher than 
the ionisation fraction predicted by the SRIO model. This 
means that for the earlier stages of the collapse, the gas in 
our model should be more strongly coupled to the magnetic 
field, which would mean there is more resistance to the col- 
lapse and the length of time spent at lower densities would 
increase, effectively introducing an almost static phase be- 
fore the initiation of the collapse. In Section [7. 1| we investi- 
gate the effect of including such a static phase, and we find 
that it only has an effect on the resulting line profiles if the 
static phase exists for > 10^ yr. Once the collapse is under- 
way, the difference in ionisation fractions between the two 
models reduces: At 5 x 10^ yr, our model predicts ionisation 
fractions of approximately 2.5 x larger than those of SRIO, 
and at 1 X 10® yr, for radii Sj 4 x 10^^ cm, the ionisation 
fractions of our model are a factor of 1.3 — 2x smaller than 
those of SRIO. For radii > 4 x 10^^ cm, the ionisation frac- 
tions of our model are factor of 2 — 4x smaller than those of 
SRIO. Since a factor of 4 difference in the ionisation fraction 
might be significant, as it means the gas in our model would 
be less strongly coupled to the magnetic field compared to 
the gas in SRlO's model and may therefore have a higher 
infall velocity, in Section |7.3| we test what happens to the 
CS line profiles when the velocity of the gas at large radii 
(> 4 X 10^^ cm) is increased. 

Note that we find that our chemical model predicts ion- 
isation fractions consistent with observations for late times: 
The estimated ionisation fractions (with respect to H nu- 
clei) of the pre-stellar cores B68 and L1544 are ~ 5 x 10~^ 
(for the central regions where Ay > 5) and ~ 5 x 10~^° (for 
radii ^ 250 0AU) respectively ( [Maret fc Bergin|2007[ |CaseUi| 
et al. 112002 1, which agree reasonably well with our estimate 



of ~ 2 X 10"" in the core centre (sC lOOOAU). 

2.2 Inside-out (collapse expansion wave) collapse 

Several solutions have been derived for the collapse of a self- 
gravitating isothermal sphere (e.g. [Larson] (|1969[); |Penston| 
( |1969[ ); |Hunter| ( |1977[ ); [Whitworth fc Summers] ( |1985| ). Here 
we concentrate on the commonly used solution derived by 
Shu (1977), which was solved by finding similarity solutions 
to the continuity (mass conservation) equation and force 
equation for an ideal isothermal fiow. The collapse expan- 
sion wave solution is a particular solution to this problem, 
which starts with the initial density distribution of a sin- 
gular isothermal sphere and with an initial velocity of zero 
everywhere. This initial configuration is in unstable equilib- 
rium, and so a small perturbation at t = yr in the core 
centre can cause the central regions to begin collapsing. The 
regions outside the collapsing centre have an unchanged den- 
sity and gravitational field, so they can maintain hydrostatic 
equilibrium and are unaffected. Material immediately above 
the collapsing inner region begins to fall, once the material 
beneath it has given way. In this way, a 'collapse expansion 
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wave' (CEW) is formed, where material inside the CEW is 
falling, and material outside is static. The CEW moves at 
the isothermal sound speed, a, which is the maximum speed 
that the collapsing inner regions can communicate with the 
static outer regions. The position of the head of the CEW is 
therefore given hy r — at at time t. The tail of the CEW is at 
the centre, where the material is approaching free-fall. This 
solution is valid for a singular isothermal sphere extending 
to infinity. For a bounded sphere, the solution breaks down 
once the CEW hits the boundary. 

As in the numerical example given by Shu (1977), we 
have used a core mass of 0.96 Mq , a core radius 1.6 x 10^^ cm 
and an isothermal sound speed of 0.2 km s~^. These values 
are applicable to Bok globules which are possible candidates 
for star formation. The collapse finishes at f = 2.5 x 10^ yr, 
when the CEW hits the boundary of the sphere. The gas 
temperature of the core, used for the chemical calculations, 
is 10 K. We ran the chemical model for a range of 'test 
shells', starting at radii of 1.0 x 10^® cm, 2.5 x 10^^ cm, 
5.0 X 10^" cm, 7.4 x 10^** cm and 1.6 x lO^'^ cm. In each 
case we calculated the density, velocity and position of the 
shell as functions of time using a spline interpolation of the 
CEW solution given in Table 2 of Shu (1977). Once a shell 
gets close to the core centre, the density and velocity rapidly 
tend towards infinity and the solution is no longer valid. At 
this point, we stop the collapse and set the velocity to zero. 
The radius, density and velocity of each shell as a function 
of time is shown in Figure [2] The point when the collapse 
of each shell halts is clearly seen by the sudden drop to 
zero in velocity. This point occurs once each shell reaches 
a distance of less than ~ 1 x lO'^^ cm from the centre of 
the core. The outer shell, at 1.6 x 10^^ cm remains static 
and at a constant density throughout the collapse. This is 
because the solution ceases to be valid once the CEW has 
reached the outermost shell, so we stop the entire collapse 
at this point, which in our model occurs at ~ 2.5 x 10^ yr. 
Table [T] shows the initial radius and density for each test 
shell for which we derive the chemical evolution, for both 
the ambipolar diffusion and inside-out collapse models. 

Note that the core mass and radius for the ambipolar 
difltusion model (A'/ = 6.80 Mq and R = 8.45 x lO^'^ cm) 
are much larger than that of the inside-out collapse model 
(Af = 0.96 Mq and R=1.6x 10^'' cm). This is because the 
ambipolar diffusion model is starting from a much earlier, 
more diffuse state, and a smaller fraction of the enclosed 
mass in the cloud would be expected to go into the protostar. 
In fact, the ambipolar diffusion model is evolving towards 
the initial state of the CEW model, and once a sufficiently 
centrally condensed core has formed it will begin collapsing 
in a more dynamic manner. Figure [7] shows that this seems 
to happen by ~ lO'^ yr, when there is a sharp increase in 
density and velocity for radii < 10^^ cm. It should also be 
noted that the collapse timescale is approximately one order 
of magnitude faster for the inside-out model than for the 
ambipolar diffusion model. 



3 CHEMICAL MODEL 

For each of the collapse models, we follow the chemical evo- 
lution of several 'test' shells spanning the radius of the core, 
in order to gain an estimate of the abundance profiles of var- 
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Figure 2. Physical evolution (radial position, velocity and den- 
sity) of the test shells in the inside-out collapse model. The grey 
dotted line, grey solid line, black dashed line, black dotted line and 
black solid line refer to the shells which start at Rq = 1.6 X 10^^, 
7.4 X IQiS, 5,0 X IQl'^, 2.5 x IQl^ and 1.0 X IQiS cm respectively. 



ious molecules at different times. The chemical reaction net- 
work is taken from the UMIST RATE06 database (Woodall 
et al. 2007), but is is limited to the elements H, C, N, O, S, 
He and Na. The species list consists of 127 gas-phase species 
and 40 mantle species and the reaction network contains 
1816 gas-phase reactions. The elemental abundances rela- 
tive to hydrogen nuclei were taken from Rawlings & Yates 
(2001) and are given in Table p] We assume that the ini- 
tial conditions are atomic, apart from hydrogen, for which 
we assume 90% is in the form of H2. In Section [TT] we also 
explore the effect of starting from a more evolved chemi- 
cal composition in each collapse model by holding the cores 
static for a period of time, ranging from 10'' — 10* yr, before 
allowing them to collapse. 

The chemical network of gas-phase chemical reactions 
includes two-body reactions, cosmic ray ionisation and cos- 
mic ray induced photoreactions. The cosmic ray ionisation 
rate used is the commonly used value of 1.30 x 10~^^ s^^ 

Shematovich et al. 2003 Garrod et al.|20()8 1 Photoreac- 



(e-g- 

tions resulting from the external radiation field are included 

in the reaction network, although we assume that the visual 
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Table 2. Abundances of each element relative to hydrogen nuclei 



H 


1.00 




C 


1.88 X 10" 


-4 


N 


1.15 X 10- 


-4 


O 


6.74 X 10" 


-4 


s 


1.62 X 10^ 


~7 


He 


0.1 




Na 


3.50 X 10- 


-7 



extinction is high (> 10 magnitudes), so these reactions have 
minimal importance. 

H2 formation on dust grains is included, the rate of 
which is given by the rate that hydrogen atoms stick and 
react to the grain surface (multiplied by 0.5 to take into 
account the fact that it takes two hydrogen atoms to form 
an H2 molecule): 



i?H2 = ^^\/3fcT/mH5'(H)n(H) 



(1) 



- X 2.35 X lQ'"T^''^nnS{K)n{Y{) cm"^ s"^ (2) 



where A is the grain surface area per unit volume, T is 
the gas temperature, n(H) is the number density of atomic 
hydrogen in cm~'^, and S is the sticking probability of atomic 
hydrogen, given by (T/102.0 + 1)"^ ^Buch fc Zhang„1991) . 
The value of A is very uncertain; we have chosen a value of 

1.5 X 10~^'^nH cm~^ (where nn is the total number density 
of hydrogen nuclei in cm~^), but values ca n be as low as 

2.6 X 10~^^nH cm~^ (JRawlings et al."l992i), or as high as 



2.1 X 10"^^nH cm-i ( |Duley fc WiUiams,1984| ). 



Freeze-out of species onto grains is included, proceed- 
ing at the same rate as given in Roberts et al. (2007). We 
have also included desorption of species owing to cosmic ray 
heating of grains at the rate given by Roberts et al. (2007), 
with the efficiency parameter cj) set to 10®. The value of (j) 
is very uncertain, but using a value of lO"" results in the 
number of molecules desorbed per cosmic ray impact being 
consistent with the theoretical estimate of Leger et al. (1985) 
(see Roberts et al. 2007 for more detail). We also note that 
other desorption mechanisms, such as desorption by H2 for- 
mation on grains and desorption by the cosmic ray induced 
UV field may also operate in dark cores, but Roberts et al. 
(2007) have shown that these processes have the same qual- 
itative effect by simply varying the value of the efficiency 
parameters, hence it is not necessary to include them all in 
this model. 

Apart from H2 formation, the only surface reactions 
that occur are hydrogenation of certain unsaturated species, 
and dissociative recombination of molecular ions. These re- 
actions are assumed to occur instantly and the products 
remain on the grain surface until they are desorbed. 

The temperature is assumed to remain at 10 K through- 
out the collapse, and the density as a function of time for 
each shell is calculated by the dynamical models described 
in Sections mi] and [O 



4 CHEMICAL RESULTS 

4.1 Ambipolar diffusion model 

Figure^ shows the fractional abundances of CS and HCO"'' 
for each test shell as they evolve with time, for the ambipolar 
diffusion model. These molecules were chosen because they 
are commonly used to search for the blue asymmetry. Ideally, 
we would have liked to model H2CO as well, which is another 
common infall tracer. However, the main formation route of 
H2CO is believed to be via the hydrogenation of CO on 
grain surfaces and subsequent desorption into the gas-phase 



(Fuchs et al. 2009). In our chemical model we have only 
included a very basic surface chemistry, which yields very 
low H2CO fractional abundances (< a few times 10"^"). In 
order to obtain a more realistic H2CO abundance, sophis- 
ticated surface chemistry would have to be implemented in 



the model, such as that employed in Cuppen et al. ( 2009 1. In 



that paper they used a Monte-Carlo approach and included 
layering of the ice, and found that the results were signifi- 
cantly different to the results obtained using a simple rate 
equation approach. This is something that we can address 
in future work, and we therefore leave modelling the H2CO 
line profiles until we have a suitable model. 

For both CS and HCO^, the chemical evolution in each 
shell is the same up to ~ 2 x 10® yrs. This is because, as 
shown in Figure fT] significant density increases only begin 
after this time. 

Up until ~ 5 X 10® yrs, there is a general increase in 
the abundances as molecules form from the atomic (with 
the exception of molecular hydrogen) initial conditions. Af- 
ter 5 X 10® — 1 X 10'' yrs, the abundances start decreasing 
due to depletion. However, including freeze-out in the model 
actually allows HCO^ to reach a peak abundance a factor 
of 5 times higher compared to when the model is run with- 
out freeze-out, before it eventually depletes. This sort of 
behaviour has been noted in previous studies (eg. Rawlings 
et al. 1992) and is complicated. It can partly be explained 
by considering the destruction channels of HCO'*': 



e + HCO+ 

H2O + HCO+ 

C -f HCO+ 



CO + H 

CO + HgO^ 

CO-I-CH+. 



(3) 
(4) 
(5) 



When freeze-out is included, the abundances of H2O and C 
decrease as they freeze onto grains. The electron abundance 
also decreases due to the depletion of easily ionised species 
such as C and Mg. Therefore the destruction rate of HCO"'' is 
initially reduced by freeze-out, allowing it to reach a higher 
peak abundance. 

Figure H] shows the radial distribution of the fractional 
abundances of CS and HCO^ across the core at 10® yr (solid 
line), 5 x 10® yr (dotted line) and 10® yr (dashed line), for the 
ambipolar diffusion model. These abundance profiles were 
generated by taking the radius of each shell at a particular 
time from Figure [l] and taking the corresponding chemical 
abundance of the shells from FigurelS] The plot clearly shows 
the contraction of the core, with the shells moving inwards 
as time progresses. 

At 1 X 10® yrs and 5 x 10® yrs, the abundances are rel- 
atively uniform across the core. At 1 x IQp yrs, when the 
collapse is well underway and there is a large density gra- 
dient across the core, both molecules suffer severe depletion 



6 J. F. Roberts, J. M. C. Rawlings and H. A. Stace 



C/3 

u 

O 



-7.5 



-8.0- 



-9.0- 



-9.5 



-10.0 

4 







- 
- 




yp 


t\ 


- 


/ 


'f - 


-. [y 










i 


- 


111 


4 

T:- 



5.0 5.5 

log[t/yrs] 



6.0 



O 



-8.0 



-8.5 



-9.0 



O 

^ -9.5 



-10.0 

-10.5 

-11.0 

4 




5.0 5.5 

log[t/yrs] 



6.0 



Figure 3. Fractional abundances of CS and HCO^ in each of the shells as a function of time for the ambipolar diffusion model. The 
black solid line, black dotted line, black dashed line, grey solid line, grey dotted line and grey dashed line are for the shells beginning at 
Ro = 8.45 X 10^'^, 7.73 x lO^'^, 6.77 x lO^'', 5.41 x 10^'', 3.93 x 10^'^ and 1.82 x 10^'^ cm respectively. The vertical dashed lines mark the 
times for which we generate line profiles. 



^ -10 

u 



o 



-12 



-14 



^ -^■■■•°-:%-€>-* 



X- 



^■ 



-¥•- 



-fl X ) l ( 



2 4 6! 

Radius/1. OxlO"cm 



10 





-8 












-9 




(^ 

1 


-X- 


■□■■•■X---0--)K : 




-10 




! 






, — , 


1 '^ 






o 
u 


-11 


1 


1 

t 
1 




^ 




-12 


/ 
-/ 






: 




-13 




,1 




-_ 



2 4 6 8 10 

Radius/ 1.0x10 "cm 



Figure 4. Fractional abundances of CS and HCO"*" of as a function of radius for the ambipolar diffusion model. The solid line is for 10^ 
yr, the dotted line is for 5 X 10^ yr, and the dashed line is for 10® yr. The asterisks, diamonds, crosses, squares, pluses, and triangles 
refer to the shells which start at Rq = 8.45 X IQi'^, 7.73 X IQl'^, 6.77 X lO^'^, 5.41 X lO", 3.93 x 10^^ and 1.82 x lO^'^ cm respectively 



for radii < 2 x 10^^ cm. In tlie outer parts of tlie core, tiie 
abundance of HCO^ actually increases as time progresses, 
due to tlie depletion effects described above. 



4.2 Inside-out collapse 

Figure Isl sliows tiie fractional abundances of CS and HCO"'" 
for eacii test sliell as tiiey evolve witli time, for the inside- 
out collapse model. The molecules show much more compli- 



cated behaviours than for the ambipolar diffusion model for 
a number of reasons, such as: (i) the density of each shell 
actually experiences an initial decrease (see Figure[2|, rather 
than monotonically increasing as in the ambipolar diffusion 
model; and (ii) the initial density distribution is a singu- 
lar isothermal sphere (rather than a uniform sphere in the 
ambipolar diffusion model), so the initial conditions of the 
shells cover a much larger range of densities. 

For HCO''' the abundance in the static outer shell (plot- 
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ted in the black solid line) sliows a very difTerent evolution 
to tlie dynamically evolving inner shells. For the inner shells, 
the molecules exhibit a large increase in abundance before 
they freeze out, due to the rapid density increase. Obvi- 
ously, for outer shells, this increase in abundance occurs 
at later times. The innermost shell collapses very quickly 
and reaches its final density of ~ 1.1 x lO'^ cm~^ after only 
~ 3.3 X lO* yr. Equilibrium between freeze-out and des- 



orption processes rapidly follows, resulting in extremely de- 
pleted abundances. 

The inclusion of freeze-out greatly increases the peak 
abundance of HCO^ before it depletes, compared to when 
the model is run without freeze-out (see Section 4.1 1, as was 
found for the ambipolar diffusion model. 

Figure [6] shows the radial abundance distribution of the 
same species at IQ"*, 5x 10*, 10^ and 2.5 x 10^ yr. For HCO+, 
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the abundances at radii > 0.5 x 10^^ cm increase with time. 
This is due to a combination of atomic to molecular conver- 
sion and depletion enhancement effects. Compared to the 
ambipolar diffusion model, the effects of depletion can be 
seen much earlier (~ 10'' yrs in the centre of the core), be- 
cause the dynamical timescale is much shorter. 

4.3 Comparison to existing chemical models of 
collapsing cores 

It is difficult to compare the abundances predicted by our 
chemical model to previously published chemical models of 
infalling cores, because of the varied range of parameters 
used (e.g. reaction data sets, elemental composition, density 
profiles, temperature profiles, etc.), but in this section we 
attempt to make a comparison with two reasonably similar 



studies. Shematovich et al. ( 2003 1 presented the abundance 



profiles of a coupled dynamical and chemical model of the 
collapse of starless cores in magnetised molecular clouds, and 
found that CS and HCO^ behaved in a similar manner to 
that which we have found in our ambipolar diffusion model; 
in Figures 2 and 3 of their paper they show that at the 
time when their model reaches a central density of nu ~ 
2 X 10*^ cm"^ (n(H2) = 10^ cm"^), for their 'low metal' 
initial abundance case (which was found to give the best fit 
to observations of L1544), the fractional abundances (with 
respect to hydrogen nuclei) of CS and HCO^ were ~ 3x 10~* 
and ~ 3 X 10~^ respectively, which are similar (within a 
factor of ~ 4) to the abundances we obtain in the outer parts 
of the core at a collapse time of 10® yr (whe n the central 

2 X 10^ cm 



density of th e SRIO model is - 



et al. 



Shematovich 
( 20031) also find at this time that CS and HCO+ are 



depleted in the core centre, in agreement with our model. 



Aikawa et al. ( 2003) modelled the molecular abundances 



in a core collapsing according to the Larson-Penston flow 



(Larson 1969 Penston 19691. Although they modelled the 



prestellar phase, starting from a uniform density, the core 
collapses to a central density of 3 x 10^ cm~^ on a timescale 
of 2 X 10^ yr, similar to the inside-out collapse model we 
have used. As one would expect, in our model we fi nd that 
CS and HCO^ deplete at an earlier time than in the 



et al. 



Aikawa 



(20031 model; by 10 yr their model shows no signs 



of depletion of CS or HCO , in contrast to our model at 
this time. By 2 x 10^ yr, however, the molecular profiles of 
CS and HCO^ show similar behaviour in both models, with 
abundances in the outer parts of the core of ~ 10~* and 
~ 10~^ respectively, and both molecules are depleted in the 
core centre. 



The spherical radiative transfer code SMMOL ('Spheri- 
cal Multi Mol') was used to generate the line profiles (Rawl- 
ings & Yates 2001). This is an approximate lambda-iteration 
code, which can be used for radiatively coupled clouds (as 
is likely to be the case for the early stages of core collapse), 
where the standard large velocity gradient (LVG) approx- 
imation is not appropriate. SMMOL was benchmarked for 
accuracy by [van Zadelhoff et al. ( 2002 1. At each radial point 
in the core, SMMOL calculates the level populations, the 
line source functions and the emergent line profile. SMMOL 
requires input tables for each molecule giving the density, 
fractional abundance, gas temperature, dust temperature, 
infall velocity and turbulent velocity at different radii for 
a particular time. To generate these tables, we took 100 
evenly spaced values of the radius between the innermost 
and outermost shells, and to obtain the abundances at each 
radii we interpolated between the data shown in Figures |4] 
and [6] We tested that this was sufficient resolution by run- 
ning the model with various grid sizes, and found that the 
spectra were invariant for > 70 radial grid points. Following 
Rawlings & Yates (2001), we used a convergence criterion 
of Arii/ni ^ lO"'' (where n; is the population of level i). As 
a boundary condition SMMOL assumes that the material 
within the inner radius is static and at a constant abun- 
dance. For the densities that we are considering, we make 
the reasonable assumption that the gas and dust are ther- 
mally coupled and have equal temperatureaM We have used 
a constant temperature of 10 K throughout the core, which 
has been found to be a valid approximation for pre-stellar 
cores such as L1498 and L1517B (Tafalla et al. 2004), al- 
though we also investigate other temperature profiles. For 
the background radiation field, we have used the standard 
interstellar radiation field of Mathis et al. (19831, supple- 



mented by a black body at 2.7 K for the cosmic background 
radiation. We have used a constant turbulent ve locity of 
0.12 km s~', which was the 1/e velocity width that 



Zhou et 



al. (19901 and Zhou et al. (19931 found to give the correct 
velocity width for observations towards B335. The veloc- 
ity and density profiles for the ambipolar diffusion and the 
inside-out collapse models are shown in Figures [T] and [S] 

The cores were assumed to be at a distance of 140 pc, 
the distance of the Taurus Molecular Cloud. For the tele- 
scope beam we have used the characteristics of the TRAM 
30m, centred on the source. These are listed in Table [3] At 
this distance the cores have angular sizes of ~ 800" and 
~ 150" for the ambipolar diffusion model and inside-out 
model respectively, so are ~ 30 times and ~ 5 times larger 
than the largest beam size (which is 27.5" at 89 GHz). 



5 GENERATING THE LINE PROFILES 

To convert these theoretical abundance profiles into observ- 
able diagnostics we calculate the line profiles of CS and 
HCO+. 

The transitions we consider are CS J = 5 — >■ 4 at 
244.94 GHz, CS J = 3 -> 2 at 146.97 GHz, CS J = 2 -^ 1 
at 97.98 GHz, and HCO+ J = 3 ^ 2 at 267.56 GHz, and 
HCO+ J = 1 -> at 89.19 GHz. These transitions were 
chosen because they are usually strong and are good tracers 
of the kinematics in protostellar infall sources (Zhou et al. 
1993; Evans et al. 2005). 



^ This assumption may be inaccurate for the early stages of 
the ambipolar diffusion model, where the density is < 10* cm~'^ 
throughout the core. However, at dust temperatures of 10-20 K, 
the dust emission only dominates the background continuum for 
frequencies > 100 GHz so the lower frequency lines should be 
unaffected by this assumption. 



Investigating the effects of chemistry on molecular line profiles of infalling low mass cores 9 



0.030 




^. .-^ 

-1.0 -0.5 0.0 0.5 1 
v/kms"' 



U.5 


CS(3-2) 




0.4 


'\ 


: 


0.3 


' 1 
' -,1 


: 


0.2 


i \ 


; 


0.1 


] 

i 


I- 


0.0 


/ 


V. 







1.4 


- CS(2-1) 




1.2 


', 




1.0 


;..' 




i^ 0.8 
0.6 




■.i '<: 




04 








0.2 


_ 




■ 


0.0 


.' ^^^ 


\ 


, , , , 1 , , , , 1 , , , 


, 1 , , , , 



-1.0 -0.5 0.0 0.5 l.( 
v/kms' 



-1.0 -0.5 0.0 0.5 1.0 
v/kms"' 



0.5 t 

HCO*(3-2) 

0.4 



0.3 






0.2- 



0.1 



-1.0 -0.5 0.0 0.5 1.0 
v/kms"' 



*<2 



HCO+(1-0) 



/ 



I / I 



-1.0 -0.5 0.0 0.5 1.0 
v/kms"' 



Figure 9. Line profiles of CS for transitions J = 5 — > 4, 3 — ^ 2, and 2—^1, HCO"*" J = 3 — > 2 and 1 — > for the ambipolar diffusion 
model at 1 X 10^ yr (solid line), 5 X 10^ yr (dotted line) and 1 X lO'' yr (dashed line). 



Table 3. Characteristics of the IRAM 30m telescope at the fre- 
quencies of the lines we are generating, taken from the IRAM web- 
site l ihttp://www.iram.fr/IRAMES/index.ht m>. HPBW is the 
half-power beam width and rj^ is the aperture efficiency. 



Frequency (GHz) HPBW (arcseconds) 



jya 



27.5 


0.61 


25 


0.60 


17.5 


0.55 


16.5 


0.55 


11 


0.42 


10 


0.39 


9 


0.36 



89 
98 
141 
147 
226 
245 
268 



6 RESULTS AND DISCUSSION 

6.1 Ambipolar diffusion line profiles 

The line profiles for CS J = 5 -^ 4, 3 ^ 2, 2 ^ 1 and HCO+ 
J — 3 ~> 2 and 1 — >■ 0, are shown in Figure [9] for collapse 
times of 1 X 10^ yr, 5 x 10^ yr and 1 x lO'' yr. 



In this model the thermal velocity is iith = \/kT/m{i) ^ 
0.043 km s"^ for CS at 10 K. At 1 x 10^ yr the maximum 
infall velocity is 0.038 km s~^, meaning that the line-width 
is dominated by turbulence and the core is essentially radia- 
tively coupled. At 5 x 10^ yr and 1 x 10® yr, however, the 
maximum infall velocity is 0.12 km s~^ and 0.31 km s~^ re- 
spectively. The systemic infall therefore begins to dominate 



the line-broadening, and at 10® yr the core centre and the 
core edge are radiatively decoupled because of the relatively 
large velocity gradient. 

Our model predicts that the CS J = 5 — >■ 4 lines are 
undetectably weak, and that the 3 — >■ 2 and 2 — >■ 1 lines only 
become detectable for times Js 5 x 10^ yr, with maximum 
intensities of ~ 0.4 and 1.0 K. This increase in intensity at 
these times is due to the increase in density of the core. The 
line profiles for CS J = 3 — >■ 2 and 2 — > 1 are very strongly 
self-absorbed and asymmetric. By 5 x 10^ yr these transi- 
tions have a noticeably asymmetric profile with a stronger 
blue-shifted peak. It can be seen that the absorption 'dip' 
is off-centre, and it is most redshifted at 5 x 10^ yr, where 
the minimum lies at ~ 0.1 km s~^, which is approximately 
the velocity of the material in the outer regions of the core 
at this time (see Figure [7|. Unlike the inside-out collapse, 
where the envelope material is static, this redshift in the ab- 
sorption is the main cause of the blue-skewed asymmetry in 
the line-profiles. Such off-centre absorption dips have been 
observed towards the class protostar, IRAM 04191, where 
they are also believed to be caused by extended inward mo- 
tions (Belloche et al. 2002). 

In the case of HCO"'", the J = 3 — >■ 2 line only becomes 
detectable at 10® yrs. The 1 — > line is detectable by 5 x 
10^ yrs, and its peak intensity strongly increases by a factor 
of ~ 7 by 10® yrs. This is due to the increase in gas density 
and in abundance of HCO''' for radii > 2 x 10^^ cm (see 
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Figure 7. Density and velocity profiles for the anibipolar diffu- 
sion model at 1 X 10^ yr (solid line), 5 X 10^ yr (dotted line) and 
1 X 10'' yr (dashed line). 



Figure 8. Density and velocity profiles for the inside-out collapse 
model at t = yr (solid line), lO** yr (dotted line), 5 X 10* yr 
(dashed line) 1 X 10^ yr (dot-dashed line) and 2.5 X 10^ yr (triple 
dot-dashed line). 



Figure ^ By 10^ yr, both the 3 — )> 2 and 1 — >• hnes are 
very strongly self absorbed with a clear blue asymmetry. 

[Pavlyuchenkov et al. (20031 also modelled CS and 
HCO^ line profiles in a spherically symmetric, coupled dy- 
namical and chemical model of a prestellar core whose evolu- 
tion is determined by ambipolar diffusion. Their model core 
was tailored to fit the stareless core L1544, and had a much 
larger mass (20 Mq) and radius (10^ AU ~ 1.5 x 10^* cm) 
than the model of SRIO. Due to the large differences core 



sizes and therefore molecular column densities, we cannot 
make a quantitative comparison to their results, but nev- 
ertheless, their resulting line profiles show notable qual- 
itative similarities to the ones we have derived. In Fig- 
ure 2 of Pavlyuchenkov et al. (20031, they present their 
HC0+(3 -^ 2) and CS(2 -^ 1) line profiles for 6.5 Myr, 
when the central density is tih ~ 2 x 10® cm"''. Comparing 
these to our Figure^ for a collapse time of 1 Myr (by which 
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Radius/10 cm 



line in the model of Pavlyuchenkov et al. ( 2003 1 has a much 



Pavlyuchenkov et al. ( 2003 1 noted that their model pre- 



in the models of Pavlyuchenkov et al. ( 2003 1 and SRIO re- 



Figure 10. The temperature profiles tested. The solid line shows 
the temperature profile of the Class source B335 (Zhou et al. 
1990), normalised to a radius of 5.9 X 10^^ cm. The dotted and 
dashed lines show the two linear temperature profiles. 



time the central density has reached ~ 10^ cm~^), we see 
that both models predict a very strong blue asymmetry in 
both lines, and a red-shifted abosrption dip. The main differ- 
ence with the results of [Pavlyuchenkov et al.| ( |2003| ) is that 
the absorption dips are somewhat wider, so the red-peaks 
therefore lie at a higher velocity (~ 0.5 km s~^) than we 
have found (0.2 - 0.3 km s~^). Furthermore, the CS(2 -)■ 1) 



higher blue-to-red peak intensity ratio; they predict a ratio 
of ~ 4, whereas we find it only to be 2.2. These differences 
appear to be because the infall velocity in the core envelope 



of the Pavlyuchenkov et al. ( 2003 1 model reaches higher ve- 
locties (~ 0.2 km s~^) compared to those in SRlO's model 
(~ 0.14 km s~^). SRlO's model reaches a higher velocity 
in the core centre (~ 0.4 km s~^), but the molecules are 
strongly depleted in this region so it does not have much 
effect on the line profiles. 



dieted too high infall velocities to fit the observations of 
L1544 by |Lee et al.| ([200T| and |Lee et al.| (|2004|). In this 



respect the model of SRIO appears to perform slightly bet- 
ter, despite the fact that the mass-to-fiux ratios are simi- 
lar in both models (~ 2 and ~ 3 times the critical value 



spectively). The differences in infall velocities could be due 
to differences in boundary conditions, and the fact that 



Pavlyuchenkov et al. ( 2003 1 do not take into account the 



drag force exerted by charged dust grains. 



6.1.1 Varying the kinetic temperature profiles for the 
ambipolar diffusion collapse 

So far, we have assumed a constant temperature of 10 K 
across the core, but now we look at how the line profile 
shapes are affected by other temperature profiles to test the 
sensitivity and possible degeneracies of the results to the 
free parameters. 

The first temperature profile we consider is that de- 



Table 4. The ratio of blue-to-red peak intensities for different 
temperature profiles for the ambipolar diffusion model at 1 X 
10^ yr. The 5 — >■ 4 line has not been included because it does not 
have an obvious double peak structure. 

Temperature profile CS(3 -s> 2) CS(2 -s> 1) 



10 K 


2.97 


2.15 


Class 


3.85 


2.17 


12K centre - 8K edge 


3.63 


2.31 


8K centre - 12K edge 


2.40 


1.99 



duced from dust continuum observations of the Class 



source B335 (Zhou et al. 19901. The shape of this profile 



comes from heating by the internal protostar in the centre, 
and heating by the interstellar radiation field at the edge 
of the core. Although the ambipolar diffusion model mostly 
describes the pre-stellar phase, by 10® yr the conditions are 
similar to the protostellar regime and therefore the class 
temperature profile may be more appropriate than the con- 
stant temperature profile at this stage. The radius has been 
normalised to 5.9 x 10^^ cm, to fit the radius of our model 
core at a collapse time of lO'' yr. Also, to understand more 
about the behaviour of the line profiles, we have tested two 
ad hoc linear temperature profiles, one of which is warm 
in the centre (12 K) and cooler on the outside (8 K), and 
the other of which is cold (8 K) in the centre and warmer 
(12 K) on the outside. The temperature profiles are shown 
in Figure [lO] 

The results are shown in Figure [ll] for the CS J — 
5 -^ 4, 3 -> 2 and 2 -)■ 1 lines at 1 x 10® yr. The overall 
characteristics of the 3 — >■ 2 and 2 — >■ 1 line profiles do 
not change (i.e. a strong blue asymmetry with an off-centre 
absorption dip) , although the ratios of the blue-to-red peak 
intensities (given in Table m are relatively sensitive to the 
temperature profiles. For the Class temperature profile and 
the linear temperature profile which is 12 K in the centre, the 
5 -> 4 line shows some absorption against the background 
continuum. The intensities are extremely low, however, so 
this feature would not be observable. 



6.2 Inside-out line profiles 

Line pro files are plotted for the same transitions as in Sec- 
tor the inside-out collapse model at 10'', 5 x 10*, 



tion 



6.1 



12 



10 and 2.5 x 10"" yr. These are shown in Figure 

For t ^ lO** yr, the inner material in the core reaches 
velocities of ~ 1.0 km s~', but the outer material remains 
static, meaning that the central regions of the core become 
largely radiatively decoupled from the outer static envelope. 
The line profiles for CS, show double peaked pro- 
files that are approximately symmetric until the final time, 
2.5 X 10^ yr. This asymmetry develops because once the 
CEW has reached the boundary of the core there is no 
longer a static outer envelope. The outer material responsi- 
ble for the absorption is now redshifted as it falls towards 
the core, and the profile kinematics are then similar to what 
was seen in the ambipolar diffusion case: the minimum in 
the absorption dip is slightly redshifted at this time (only), 
at ~ 0.02 km s~'. Thus, it would appear that if there is 
no imposed kinetic temperature gradient then, regardless of 
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Figure 11. Plot of CS spectra at 1 X 10^ yr for different temperature profiles, for ambipolar diffusion collapse. Top row: constant 
temperature of lOK; second row: Class temperature profile; third row: linear temperature profile which is 12 K in the centre and 8 K 
at the edge; bottom row: linear temperature profile which is 8 K in the centre and 12 K at the edge. 



the nature of the collapse, the presence of a blue asymmetry 
is indicative of an infalling outer envelope. 

The peak intensities of all three CS lines decrease with 
time. Referring to FigurelGl this appears to be a consequence 
of the decrease in CS abundance in the central regions of the 
core due to freeze-out, as well as the decrease in density (at 
a fixed radius) which is shown in Figure [s] which affects the 
collisional excitation as well as decreasing the CS column 
density. 

For early collapse times (^ 1 x 10^ yr), the CS line in- 
tensities for the inside-out collapse model are much higher 
than for the ambipolar diffusion model, because at this stage 
the inside-out collapse model has much higher densities and 
there is still a relatively high CS abundance in the core 
centre. By 2 x 10^ yr in the inside-out collapse model and 
1 X 10® yr in the ambipolar diffusion model, the central core 
densities and CS depletion are similar in the two models, 
which is reflected in the similar CS line intensities at these 
times. 

The HCO"'' line profiles are strong, symmetric and dou- 



ble peaked. A blue asymmetry develops at later times; by 
10^ yr, the 3 — >■ 2 line has a marked asymmetry and en- 
hanced line width. At 2.5 x 10^ yr both lines are blue- 
asymmetric. Whilst the 1 — >■ line line strength is reason- 
ably constant (except at late times) , that of the 3 — >■ 2 line 
declines steadily with time. Again, as a result of the higher 
densities, the HCO'*' line intensities are much greater in the 
inside-out collapse model compared to the ambipolar diffu- 
sion model (apart from when the ambipolar diffusion model 
reaches 10® yr, when it also has a high central density). 



6.2.1 Varying the kinetic temperature profiles for the 
inside-out collapse 

In the context of the inside-out collapse model (with a static 
outer envelope) it is usually stated that, in order to have 
asymmetric line profiles with a stronger blue-shifted peak, 
a negative excitation temperature gradient is needed (Zhou 
et al. 1993). 

In Figure |13[ the thick solid line plots the excitation 
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Figure 12. Line profiles of CS for transitions J = 5 — ^ 4, 3 — > 2, and 2—^1, HCO"^ J = 3 —^ 2 and 1—^0 for the inside-out collapse 
model at 1 X 10* yr (solid line), 5 X 10* yr (dotted line), 1 X 10^ yr (dashed hne) and 2.5 X 10^ yr (dot-dashed line). 
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Figure 13. The excitation temperature, Tex, of the CS(3 — > 2) 
line as a function of radius, for the inside-out collapse model at 
a collapse time of 10^ yr. The thick solid line shows Tex for a 
constant kinetic temperature profile of 10 K, using the CS abun- 
dance profile of Figure[6] The thin lines show Tex using a constant 
CS abundance of 2.5 X 10~*, with a constant kinetic temperature 
gradient of 10 K (thin solid line), a Class 0-type temperature pro- 
file (dotted line), a linear temperature profile with 12 K in the 
core centre and 8 K at the core edge (dashed line), and a linear 
temperature profile with 8 K in the core centre and 12 K at the 
core edge (dot-dashed line). 



temperature, Tex, of t he CS(3 — >■ 2) line for the spec- 
trum plotted in Figure 12 for a collapse time of 10^ yr. At 
this time, there is a negative excitation gradient for radii 
^ 5.2 X 10^® cm, increasing from 2.2 K at the outer radius 
to 4.8 K at 5.2 x lO^*' cm. For radii of 1.2 x 10^*^ < R < 
5.2 x 10^^ cm, however, there is a positive excitation gradi- 
ent. The infall radius at this time is ~ 6 x 10^® cm, so the 
infalling material is within the region with the positive exci- 
tation temperature gradient. Therefore, the central infalling 
blue-shifted material does not have an excitation temper- 
ature higher than the red-shifted infalling material at the 
edge of the collapse expansion wave, and this can explain 
why we do not see the blue asymmetry at this time. We 
note that for radii < 1.2 x 10^^ cm, the excitation tempera- 
ture sharply increases towards the core centre, but for such 
small radii the CS is heavily depleted. 

If a negative kinetic temperature gradient is imposed, 
this can change the excitation temperature gradient in the 
core, and could give rise to the blue asymmetry. However, 
we also found that the CS line profiles remained symmetric 
for all times SC 1 x 10^ yr even when adopting the various 
temperature profiles described in Section [6.1.1| Rawlings & 
Yates (2001) also found that CS exhibited symmetric line 
profiles for inside-out collapse, and attributed this to the 
depletion of CS in the central regions. 

To test this theory, we have calculated the line profiles 
for CS at 10^ yr using the temperature profiles described in 
Section [6. 1.1[ but this time assuming a constant fractional 
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Table 5. The ratio of blue-to-red peak intensities for different 
temperature profiles for the inside-out collapse at 1.0 X 10^ yr, 
assuming a constant CS abundance across the core. 



Temp, profile 



CS(5 -5> 4) CS(3 -s> 2) CS(2 -s> 1) 



10 K 1.29 

Class 1.52 

12K centre - 8K edge 1.29 

8K centre - 12K edge 1.10 



1.03 
1.10 
1.08 
0.95 



1.01 
1.03 
1.04 
0.98 



abundance of CS across the core. We have used a fractional 
abundance for CS of 2.5 x 10~*, which is the CS abundance 
calculated for the outer radius of the core in the original 
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model at 10 yr. The line profiles are plotted in Figure 
and the ratios of blue-to-red peak intensities for each tem- 
perature profile are given in Table [5] The excitation temper- 
atures of the 3 — > 2 transition are plotted in Figure |13[ 

From Figure[l4]we can see that all of the CS transitions 
now show high velocity wings since CS is not depleted in the 
high velocity central core regions. For all of the temperature 
profiles, the 2 — > 1 lines are more or less symmetric, but the 
5 — >■ 4 and 3 — >■ 2 lines (which trace the denser, faster mov- 
ing inner material) show stronger asymmetries. The red-blue 
asymmetries in the lines (although not strong) broadly re- 
flect the expected dependency on the excitation temperature 
profile. The class temperature profile gives the strongest 
blue asymmetries, which can be explained by the very high 
excitation temperature in the core centre (where there is 
still a high abundance of CS since we have neglected deple- 
tion). It is also interesting that even for the constant kinetic 
temperature profile, the 5 — >■ 4 line shows a blue asymmetry. 

These results support the interpretation of Zhou et al. 
(1993), that a kinetic temperature gradient with a warm 
core centre can explain the blue asymmetry for inside-out 
collapse. However, they also indicate that, in order for CS 
to show the blue asymmetry, it cannot be depleted in the 
core centre. 



7 DISCUSSION 

It has been demonstrated in the previous section that the 
abundance profile of the molecules across the core (in partic- 
ular whether or not there is depletion) can have a significant 
impact upon the line profiles. We therefore first discuss the 
sensitivity of the models to different chemical conditions. 



7.1 Initial molecular abundances 

We have assumed that, in each of our models, the core be- 
gins its collapse from atomic initial conditions (apart from 
hydrogen which is initially 90% molecular). This, however, 
might not be accurate, especially for the inside-out collapse 
model which starts from a more evolved state. This has been 
recognised in previous studies. For example, Rawlings et al. 
(1992) assume that the core is initially at a uniform density 
of 2.8 X 10'^ cm~^, with the initial abundances taken from 



a previous model by Nejad et al. (19901, and then allow 
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Figure 15. Line profiles of CS for transitions J = 5 — > 4, 3 — > 2, 
and 2 — > 1 for the ambipolar diffusion model at 1 X 10^ yr (solid 
line), 5 X 10^ yr (dotted line) and 1 X 10^ yr (dashed line). In this 
model the core has been static for tstat = 10^ yr before it begins 
collapsing. 
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modified free-fall collapse to the isothermal sphere config- 
uration, where the inside-out collapse eventually begins. In 



Figure 16. Line profiles of HCO+ for transitions J = 3 — > 2 
and 1 — > for the ambipolar diffusion model at 1 X 10^ yr (solid 
line), 5 X 10^ yr (dotted line) and 1 X lO'' yr (dashed line). In this 
model the core has been static for tstat = 10^ yr before it begins 
collapsing. 



Rawlings & Yates (2001), they also consider the possibility 
that once the core has reached the isothermal sphere state, 
it can remain in pressure-balanced equilibrium prior to the 
inside-out collapse. 

It is debatable what the initial abundances for such col- 
lapse models should be, since the pre-stellar core lifetime 
is not yet well determined ( [Andre et al.|[2009| |. We have 
therefore investigated the effect on the line profiles for both 
models of keeping the core static for a period of time, tstat, 
ranging from 10^ to 10^ yr, allowing the chemistry to reach a 
more evolved state before the collapse is initiated. The draw- 
back of this investigation for the inside-out collapse model is 
that it may be unrealistic that a core would be static in the 
isothermal sphere configuration preceding the onset of the 
CEW. Instead, it is more likely that it is contracting from 
lower densities, and therefore the chemical timescales could 
be longer than those presented here. 

For the ambipolar diffusion model, we find that there 
is no effect on the line profiles presented in Figure [9] for 
tstat ^ 10* yr. For tstat = 10^ yr, there is an effect on the 
CS line profiles, shown in Figure [Ts] In Figure |9] the peak 
intensities of the CS(3 — >■ 2) and (2 — > 1) lines increase as the 
core evolves from 5 x 10^ yr to 10^ yr, but for tstat = 10^ yr 
we find that those lines reach their maximum peak intensity 
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Figure 14. Plot of CS spectra for 1.0 X 10® yr for different temperature profiles, for inside-out collapse assuming a constant CS abundance 
across the core. Top row: constant temperature (lOK); second row: Class temperature profile; third row: linear temperature profile 
which is 12 K in the centre and 8 K at the edge; bottom row: linear temperature profile which is 8 K in the centre and 12 K at the edge. 



at 5 X 10® yr, and then they decrease by 10® yr. This is 
because for i^tat = 0, CS reaches its peak abundance between 
5.6 X 10® yr (for the inner sheU) and 7.0 x 10® yr (for the 
outer shell), but for fstat ~ 10® yr, the peak abundance is 
reached at slightly earlier times of 5.2 x 10® to 6.0 x 10® yr. 
Therefore for iatat = 10® yr, by a collapse time of 10® yr, the 
CS abundance in all shells is already past its peak, resulting 
in a lower abundance across the core compared to tgtat = 0. 
For fstat 5S 10® yr, CS is severely depleted, so the peak 



> 10 



7 



yr, 



CS 



intensities are very much reduced. For ts 
becomes undetectable. 

For HCO"^, the line profiles are not affected until istat ~ 
10® yr. For this iatat, the initial HCO^ abundance for the 
collapse stage is very high (3.2 x 10"^ compared to ~ 10"^" 
for tstat = 0), resulting in much higher peak intensities for 
the earlier collapse times (10® and 5 x 10® yr), compared to 
those with fstat = 0. 

For fstat ~ W^ yr, shown in Figure 
dance of HCO^ is even higher (8.7 x 10~ 
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the initial abun- 
), and line inten- 



By fstat = 10 yr, however, HCO^ begins to deplete be- 
fore the collapse begins (the initial HCO"*" abundance for the 
collapse stage is now only 1.3 x 10^') , and the line intensities 
decrease compared to fstat = 10^ yr. 

For the inside-out collapse model, there is no effect on 
the CS or HCO+ hue profiles until fstat = 10® yr. For this 
fstat, even for a collapse time of 10* yr, the CS becomes 
severely depleted in the inner shells (for the three innermost 
shells the abundances of CS are 4.9 x 10"^®, 5.6 x 10"^^ 



and 1.2 x 10"^ for fs 



10 yr, compared to 7.6 x 10 



1.2 X 10~** and 1.3 x 10~* for fstat = 0). This results in large 
decrease in the line intensity during the collapse for all the 
CS lines considered. This effect is particularly severe for the 
CS(5 — >■ 4) which traces the innermost parts of the shell; for 
a collapse time of lO'' yr the line intensity reduces from 1.2 
for fstat = to 0.12 for fstat = 10® yr. 

For HCO+, for fstat = 10® (plotted in Figure [l7|, at all 



sities increase even further compared to fstat ~ 0. 



times during the collapse the abundances in the outer parts 
of the core are greater than the abundance in the same shells 
for fstat = 0. However, at this fstat, the inner shells are more 
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Figure 17. Line profiles of HCO+ for transitions J = 3 — > 2 
and 1—^0 for the inside-out model at 1 X lO** yr (solid line), 
5 X 10* yr (dotted line), 1 X 10^ yr (dashed line) and 2.5 X 10^ 
yr (dot-dashed line). In this model the core has been static for 
tstat = 10^ yr before it begins collapsing. 



depleted. Because of this, for fstat = 10^ yr, the HCO^(3 — > 
2) lines (which trace the inner parts of the core) have lower 
peak intensities than for iatat = 0, and the HCO^(l — ^ 0) 
lines (which trace the lower density outer material), show 
higher peak intensities for tstat = 10^ yr. 

By tatat ~ 10^ yr, the CS lines have practically disap- 
peared because of depletion, and even the HCO'''(l — > 0) 
lines start decreasing in intensity as the HCO"'' depletion 
starts to affect larger radii. 

As with the ambipolar diffusion model, keeping the core 
static for a period of time aflects the CS and HCO^ lines dif- 
ferently. For both models, including a static period tends to 
weaken the CS lines due to depletion, but for moderate val- 
ues of istat (~ 10® yr for the ambipolar diffusion model and 
~ 10^ yr for the inside-out collapse model), the HCO^ inten- 
sities can increase. This is because, as explained in SectionHj 
HCO^ is a 'depletion enhanced' molecule. The CS:HCC5^ 
intensity ratios could therefore be a useful indicator of the 
initial state of chemical evolution of a collapsing core. The 
inside-out collapse model results indicate that for CS to be 
detectable during the collapse phase, the core cannot have 
been static with the isothermal sphere density configuration 
for more than 10® yr. 



7.2 Visual extinction and inclusion of 
photoreactions 

In our models we have assumed that the cores are sufficiently 
embedded within a larger molecular cloud such that the vi- 
sual extinction at the edge of the envelope is > 10 mag- 
nitudes, practically suppressing all of the direct photore- 
actions. It is possible, however, that photoreactions play a 
more significant role in the chemistry, particularly at larger 
radii where the material is less shielded. This is particularly 
likely in cases where the parent molecular cloud has been dis- 
persed (eg. the pre-stellar core B68 ( Maret fc Bergiii||2007 1 
and the Class object B335 ( |Zhou et al.||1993| )). 

To test the importance of these effects we have re-run 
our chemical model for the shells at the outer edges of 
the envelopes for both collapse models. We find that, for 
Av s; 5 mag, the abundances of CS and HCO^ are signif- 



icantly reduced. For CS, this is because it is directly de- 
stroyed by photons. The higher ionisation level suppresses 
HCO"'' through dissociative recombination. 

These effects lead to reduced molecular abundances in 
the outer parts of the core, resulting in steeper gradients in 
the abundance profiles, and thus changing the line profile 
shapes. For example, in the ambipolar diffusion model, we 
suggested that the blue asymmetries in CS and HCO^ lines 
are caused by the redshifted material at the front of the 
envelope preferentially absorbing the red emission from fur- 
ther inside the cloud. If there was less material in the outer 
envelope there would be less absorption, and we therefore 
expect a less pronounced blue asymmetry. 



7.3 Effect of ionisation fractions for the ambipolar 
diffusion model 

In Section ["2.1| we mentioned that the ionisation fraction pre- 
dicted by our chemical model is ~ 4 times less than the ioni- 
sation fraction predicted by the model of SRIO at late times 
(10® yr) , for large radii (> 4 x lO^'^ cm). Since the dynamics 
of SRlO's model, which we used as an input to our chemical 
model, depend on the ionisation fraction, it is possible that, 
to be consistent with our chemical model, the outer parts of 
the core (i? > 4 x 10^^ cm) should be less strongly coupled 
to the magnetic field and could therefore have larger infall 
velocities. 

To understand what effect this would have, we have 
recalculated the CS line profiles for a collapse time of 10® yr, 
but doubled the infall velocity for 7? > 4 x 10^^ cm. We find 
that even for such high infall velocities in the outer parts 
of the core, the blue-shifted peak of the CS line profiles is 
almost completely unaffected. The red-shifted peaks, on the 
other hand, show a decrease in peak intensity, and for the 
CS(2 — >■ 1) line the red-shifted peak is shifted further to the 
right; in the original spectrum the peak intensity is 0.51 K 
at 0.26 km s"'^, but with the increased velocities in the outer 
core, the peak intensity decreases to 0.26 K at 0.31 km s~^. 

This is the expected behaviour, since the higher velocity 
red-shifted material at the front of the core (i.e. closest to 
the observer) is now more strongly radiatively coupled with 
the red-shifted high velocity material in the core centre, and 
therefore can absorb more of the red-shifted emission from 
the core centre. Note also that in this 'high-velocity' model, 
the infall velocities are similar to those of |Pavlyuchenkov et| 
al. ( 2003 1 (~ 0.2 km s~^), and the blue-to-red peak intensity 
ratios for the CS(2 — ^ 1) line predicted by both models are 
now in very good agreement with each other (~ 4). 



8 CONCLUSIONS 

In this paper we have modelled and compared the line pro- 
files of several transitions of CS and HCO"'' at various times 
for cores undergoing (a) ambipolar diffusion regulated col- 
lapse and (b) inside-out collapse. Rather than assuming a 
constant abundance for each molecule, we have coupled a 
chemical model to the dynamical models to try to generate 
realistic, time dependent abundance profiles, which are then 
used as an input to the radiative transfer code to calculate 
the line profiles. We have found four important conclusions: 
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(i) For the inside-out collapse model to exhibit the blue 
asymmetry it is necessary to suppress freeze-out and im- 
pose a negative temperature gradient on the core, whereas 
the ambipolar diffusion model exhibits the blue asymme- 
try even when there is significant freeze-out, and for all the 
kinetic temperature profiles tested. This conclusion is im- 
portant because, for example, CS has been observed to be 
significantly depleted within dense cores (e.g. Bergin et al. 
2001; Pagani et al. 2005), and yet exhibits the blue asym- 
metric line profile (e.g. Mardones et al., 1997; Lee et al., 
1999). Zhou et al (1993) found that their observations of 
CS in B335 were very well fitted by the inside-out collapse 
model, but in these studies the CS abundance was assumed 
to be constant across the core. Given the strong evidence 
for depletion of CS in such cores, it may be worth investi- 
gating in more detail the CS abundance profile in this core 
using interferometric observations, and seeing if the inside- 
out model is still able to reproduce the CS observations with 
a realistic abundance profile. Note that Belloche et al. (2002) 
also found that the inside-out collapse model could not re- 
produce the blue asymmetry observed towards the class 
protostar IRAM 04191 unless the collapse expansion wave 
radius is very large (Js 10 000 AU), but then this resulted in 
very high velocities at small radii, giving rise to excessively 
large linewidths. They estimated that CS is depleted by a 
factor of ~ 20 in the centre of this core. Since they took into 
account this depletion in their radiative transfer model, this 
could explain why their inside-out collapse model failed to 
reproduce the observations. 

The fact that a negative kinetic temperature gradient is 
required in the inside-out collapse model but not in the am- 
bipolar diffusion model is due to the different line formation 
processes: in the ambipolar diffusion model the blue asym- 
metry arises because the foreground absorbing material is 
also infalling and redshifted, but for the inside-out collapse 
model it is because the blue emission originates from the 
central regions of the core where the excitation temperature 
is relatively high, whereas the red emission originates from 
further away from the centre where there is relatively lower 
excitation temperature. 

(ii) There are marked differences between the two models 
in the line strengths, asymmetries and time-dependencies 
of the line profiles. For the ambipolar diffusion model, the 
CS(2 — >■ 1) and (3 — >■ 2) line profiles are strong and demon- 
strate a blue asymmetry (with an absorption trough that is 
not at zero velocity). The CS(5 — >■ 4) line is not detectable. 
The HC0+(3 -)> 2) line is only detectable after ~ 10*^ years, 
whilst the (1 — > 0) line strength strongly increases with time. 
For the inside-out model, all three CS lines are strong from 
the earliest stages, but decline with time. The HC0'''(1 — > 0) 
line is also strong at all times, but the strength of the (3 — > 2) 
line steadily declines. Thus, subject to caveats about possi- 
ble degeneracies, these various characteristics may be use- 
ful in determining the dynamical structure and evolutionary 
status of low-mass infall regions. 

(iii) The initial molecular abundances in collapsing cores 
affect the evolution of the line profiles. In particular, if the 
core is held static for a period of time before the collapse, 
the peak intensity ratio between the CS and HCO'*' lines can 
be strongly affected. This is a major source of uncertainty 
in the theoretical interpretation of observed line profiles. 



This work demonstrates that the inside-out collapse 



model of Shu (19771, often referred to as the 'standard' in- 
fall model (e.g. iDunham et al.||2010| |Hung et al.|2010| ), has 
difficulty reproducing the blue asymmetry in the line pro- 
files when the chemical evolution of the core is taken into 
account, due to depletion in the core centre and a lack of ex- 
tended inward motions. The ambipolar diffusion model does 
not suffer from this problem, because even when the core 
has reached an advanced, centrally condensed state, there 
are still extended inward motions. 

Before using line profiles as diagnostics of a particular 
collapse model, it is essential to investigate the abundance 
profile of the tracer molecule, in particular to see if there is 
any evidence for depletion. 
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